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1 .  INTRODUCTION 


The  print  and  fire  processing  of  thick  film  circuits  ensures  that 
there  always  will  be  some  degree  of  chemical  interaction  between  the  film 
and  the  substrate,  because  all  common  substrate  materials  are  soluble  to 
some  degree  in  the  glasses  used  in  thick  film  inks.  This  interaction  is 
primarily  responsible  for  the  development  of  adhesion  between  the  thick 
film  and  the  substrate,  but  it  also  leads  to  changes  in  the  composition 
of  the  glass  with  the  net  result  that  the  physical  properties  of  the  glass 
will  change.  These  changes  in  physical  properties  of  the  glass  will  result 
in  modified  kinetics  for  the  various  microstructure  development  processes 
and  all  electrical  properties  of  the  resistors  are  related  to  the  micro¬ 
structure.  The  goal  of  this  research  program  is  to  develop  a  sufficient 
level  of  understanding  of  the  phenomena  involved  so  that  appropriate  models 
can  be  developed.  These  models  should  lead  to  the  writing  of  specifications 
for  impurity  limits  and  additive  ranges  for  substrates,  and  to  recommenda¬ 
tions  concerning  glass  composition  and  processing  conditions. 

Previously  reported  studies  [1-5]  under  this  program  have  established 
the  magnitude  of  the  effects  resulting  from  chemical  interaction  between  a 
thick  film  resistor  and  a  ceramic  substrate,  and  have  determined  the  speci¬ 
fic  influence  on  important  properties  of  the  resistor  glass.  The  rates  of 
dissolution  of  two  substrates,  96%  Al^  (AlSiMag  614)  and  99.5%  A^O-j 
(AlSIMag  772),  in  two  lead  borosilicate  glasses  (63  w/o  PbO-25  w/o  BgO^-l 2 
w/o  Si O2  and  70  w/o  Pb0-20  w/o  B203-10  w/o  Si 02 )  were  measured  at  various 
temperatures.  The  rate  limiting  steps  for  each  substrate-glass  system  were 
determined  in  all  appropriate  temperature  ranges,  and  analytical  equations 


-2- 


were  developed  to  predict  the  substrate  recession  as  a  function  of  time  and 
temperature  for  thick  film  resistors.  The  saturation  solubility  of  AlSiMag 
614  substrates  in  63-25-12  glass  was  determined  as  a  function  of  temperature, 
and  these  results  were  combined  with  the  dissolution  rate  studies  in  order 
to  test  various  kinetic  models.  Studies  of  the  kinetics  of  initial  stage 
sintering  of  glass  particles  were  conducted  in  order  to  determine  the  ratio 
of  surface  tension  to  viscosity  for  the  two  standard  glasses,  and  the  stand¬ 
ard  glasses  with  additions  of  substrate  ingredients.  Both  the  magnitude  and 
the  activation  energies  of  this  ratio  were  found  to  be  significantly  dif¬ 
ferent  for  the  glasses,  confirming  the  extreme  sensitivity  of  this  parameter 
to  small  changes  in  glass  composition.  The  viscosity  of  the  63-25-12  glass 
was  measured  as  a  function  of  amount  of  dissolved  AlSiMag  614  substrate 
from  450°  to  850°C.  The  isothermal  viscosity  was  found  to  increase  by  a 
factor  of  20  with  10  w/o  dissolved  substrate  relative  to  the  standard  lead 
borosilicate  glass.  The  surface  tension  of  the  63-25-12  glass  was  measured 
as  a  function  of  the  amount  of  dissolved  AlSiMag  614  substrate  from  600°C 
to  900°C.  The  surface  tension  increased  with  increasing  amount  of  dissolved 
substrate  with  the  effect  being  most  pronounced  at  the  higher  temperatures. 
The  effects  of  the  minor  constituents  (4%)  in  the  AlSiMag  614  substrate 
composition  on  viscosity  and  surface  tension  of  the  glasses  was  determined 
to  be  larger  than  predictions  based  on  the  principal  of  additivity.  The 
sheet  resistance,  temperature  dependence  of  resistance,  and  the  current 
noise  were  measured  for  thick  film  resistors  as  a  function  of  the  amount 
of  substrate  dissolved  in  the  resistor  for  various  firing  temperatures  at 
constant  firing  time,  and  for  various  firing  times  at  800°C.  Large  varia¬ 
tions  in  these  three  properties  were  observed,  and  the  changes  were 


qualitatively  correlated  with  changes  in  viscosity  of  the  glass.  The  micro 
structure  development  and  charge  transport  models  used  to  correlate  the 
results  indicated  a  retardation  of  microstructure  development  in  the 
resistors  as  the  amount  of  dissolved  substrate  increased  which  led  to  chang 
ing  proportions  of  sintered  and  non- sintered  contacts  in  the  RuC^  networks 
within  the  body  of  the  resistor.  The  non-sintered  contacts  were  modeled 
using  MIM  devices,  and  the  dominant  charge  transport  mechanism  for  glass 
films  >7  nm  at  room  temperature  and  above  was  Schottky  emission,  while 
quantum  mechanical  tunneling  became  important  for  thinner  films  and/or 
lower  temperatures. 
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2.  DENSITY  OF  GLASSES 

The  dissolution  of  AISiMag  614  (96%  A^O-j)  substrate  in  the  standard 
glass  (63  w/o  PbO-25  w/o  B203-12  w/o  Si02 )  was  shown  [2]  to  proceed  through 
molecular  diffusion  and  free  convection  as  rate  limiting  steps.  Free  con¬ 
vection  could  be  caused  by  gradients  in  temperature,  surface  tension,  or 
density  due  to  the  dissolution  of  substrate  in  the  glass.  Convection  due 
to  temperature  differences  caused  by  the  reaction  can  be  considered 
negligible  because  the  experiments  were  carried  out  at  sufficiently  high 
temperatures  (i760°C)  that  any  thermal  gradients  caused  by  the  dissolution 
would  have  been  balanced  by  heat  radiation  from  the  bulk  glass.  Surface 
tension  effects  need  not  be  considered  for  these  experiments  because  the 
substrates  were  completely  immersed  in  the  glass,  which  leaves  density 
gradients  as  the  most  probable  driving  force  for  convection. 

The  density  was  determined  as  a  function  of  composition  and  tempera¬ 
ture  in  the  gold  sphere  viscometer  [4]  utilizing  Archimede's  principle: 

Glass  density  -  ^ 

x  density  of  gold  (1) 

The  density  results  for  the  standard  glass  and  four  other  compositions 
doped  with  AISiMag  614  substrate  or  alumina  are  given  in  Fig.  1.  The 
density  (p)  for  each  glass  can  be  adequately  represented  as  a  linear  func¬ 
tion  of  absolute  temperature  (T)  in  the  range  measured.  The  densities  in 
gm/cc  are: 

Standard  glass:  p  *  5.2607  -  1.1068  x  10"^  T  (2) 

6  w/o  sub.  glass:  p  *  4.6997  -  7.6200  x  10“4  T  (3) 


Temperature  (  C) 

Figure  1.  Variation  of  Glass  Density  with  Composition  and  Temperature. 
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10  w/o  sub.  glass:  p  =  4.5096  -  6.6000  x  10"4  T  (4) 

14  w/o  sub.  glass:  p  =  4.1457  -  3.8800  x  10'4  T  (5) 

8  w/o  A1203  glass:  p  =  4.7163  -  7.1400  x  10"4  T  (6) 

The  decrease  in  density  of  the  glass  due  to  substrate  dissolution  is  suf¬ 

ficiently  large  to  lead  to  density  gradients  which  could  be  responsible  for 
the  natural  convection  following  the  molecular  diffusion  process  during  the 
dissolution  of  the  substrate.  The  temperature  coefficient  of  density,  for 
the  temperature  range  investigated,  is  seen  to  decrease  with  increasing 
additions  of  substrate. 

A  comparison  of  the  density  of  the  8  w/o  A1203  glass  with  the  substrate 
glasses  in  Fig.  1  shows  that  the  minor  constituents  (3.6%  Si02,  1.3%  MgO, 

0.1%  CaO,  0.2%  Fe203)  in  the  AlSiMag  614  substrate  have  a  significant  ef¬ 
fect.  The  density  of  a  liquid  may  be  approximated  to  be  an  additive  pro¬ 
perty  of  its  constituents,  but  nonlinear  behavior  would  be  expected  if  the 
constituents  differ  greatly  in  density  or  if  there  is  a  change  in  the  struc¬ 
ture  of  the  liquid  due  to  the  addition.  Additivity  factors  for  room  tempera¬ 
ture  density  based  on  percent  by  weight  of  PbO,  B203,  Si02,  and  A1203  have 
been  reported  [6]  to  be  9.6,  1.9,  2.3,  and  4.1,  respectively.  Addition  of 
A1 203  to  a  high  lead  glass  would  be  expected  to  decrease  the  density,  and 
addition  of  Si02  would  bring  the  density  down  even  further.  Thus  the  addi¬ 
tivity  factors  qualitatively  explain  the  fact  that  the  addition  of  substrate 
lowers  the  density  more  than  the  addition  of  an  equal  amount  of  pure  A1203 
However,  the  differences  cannot  be  quantitatively  accounted  for  by  consider¬ 
ing  additivity  factors  alone;  possible  differences  in  glass  structure,  as 
indicated  by  the  anomalous  behavior  of  the  viscosity  of  8  w/o  A1 203  glass 
[3],  may  be  responsible  for  the  observed  anomaly  in  density. 
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3.  Ru02  RIPENING 

The  experimental  techniques  and  some  preliminary  results  for  studies 
of  the  rate  of  dissolution  of  small  particles  and  growth  of  large  particles 
(ripening)  of  Ru02  in  glass  were  previously  presented  [4].  The  complete 
results  of  surface  area  measurements  on  Ru02  samples  held  for  varying  times 
in  the  standard,  2  w/o,  4  w/o,  and  10  w/o  substrate  doped  glasses  at  950°C 
are  given  in  Fig.  2.  The  extensive  measurements  of  Ru02  particle  growth 
in  the  standard  and  10  w/o  substrate  glasses  showed  that  after  a  rapid 
decrease,  the  surface  area  decreased  very  slowly  at  longer  times.  The 
results  also  indicate  that  the  ripening  process  is  considerably  retarded 
by  the  addition  of  substrate  to  the  standard  glass.  The  amount  of  Ru02 
relative  to  the  glasses  was  varied  from  20  to  30  w/o  for  the  standard  and 
10  w/o  substrate  doped  glasses,  and  the  effect  of  this  variation  on  the 
ripening  kinetics  is  seen  to  be  insignificant.  For  ripening  studies  at 
888°  and  1000°C,  composites  made  of  Ru02  and  the  standard, 4  w/o,  10  w/o, 
and  14  w/o  substrate  doped  glasses  were  used  in  order  to  cover  a  wider 
glass  composition  range. 

If  diffusion  controlled  solution-precipitation  is  the  rate  limiting 
step  for  the  process  of  ripening  of  solid  particles  dispersed  in  a  continu¬ 
ous  matrix,  the  average  particle  radius  obeys  the  following  time  dependence. 

[r  (t)]3  -  [r  (o)]3  =  [f  Co^-fV°  °]  t  (7) 

In  this  equation,  r  (t)  is  the  average  particle  radius  at  time  t,  r  (o)  is 
the  average  particle  radius  time  zero,  CQ  is  the  equilibrium  solubility, 
YS1  is  the  solid-liquid  interfacial  energy,  VQ  is  the  molar  volume,  D  is 
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the  diffusion  coefficient,  and  R  is  the  gas  constant.  If  the  particle  growth 
occurs  by  a  phase  boundary  reaction  controlled  solution-precipitation  process, 
the  average  particle  radius  increases  according  to 

C?  (t)]2  -  [F  (o)]2  *  [(f)2  C°-Y|}<°  **]  t  (8) 

where  K^.  is  the  transfer  coefficient.  Lease  squares  fits,  of  the  experimental 
data,  to  these  equations  for  varying  glass  compositions  at  three  different 
temperatures  indicated  that  the  interface  reaction  is  the  slowest  step  of  mass 
transport  in  this  system.  Plots  of  ([r  (t)]  -  [r  (o)]  )  versus  time  are 
are  straight  lines  for  all  the  three  temperatures  investigated  as  shown  in 
Figs.  3-5.  The  slopes  of  the  lines,  which  are  proportional  to  C0yS]Kj,  are 
significantly  reduced  by  the  addition  of  substrate  to  the  standard  glass. 

The  values  of  Cqys1Kj  were  calculated  from  the  slopes,  and  least  squares 
analyses  on  the  data  to  an  Arrhenius  fit  were  performed  to  give  the  follow¬ 
ing  equations: 

Standard  glass:  Cqys1Kt  =  8.577  x  106  exp  [-(51601 

±  6416)/T]  (9) 

4  w/o  sub.  glass:  Cqys^Kj  3  6.533  x  10^  exp  [-(49187 

±  7559)/T]  (10) 

10  w/o  sub,  glass:  Cqys1Kt  *  9.206  x  104  exp  [-(47637 

±  6796)/T]  (11) 

14  w/o  sub.  glass:  Cqys1Kt  3  5.305  x  105  exp  [-50085/T]  (12) 

o 

In  these  equations,  CQ  is  in  gm/cc,  ys^  in  0/cm  ,  and  Ky  in  cm/s.  An  aver¬ 
age  exponential  term  of  -49626/T  was  selected  to  represent  all  four  compo¬ 
sitions  because  the  differences  among  the  individual  glasses  were  within 
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Figure  3.  Ru02  Ripening  in  Glasses  at  888°C. 
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Figure  4.  Ru02  Ripening  in  Glasses  at  950°C. 


Rigure  5.  Ru02  Ripening  in  Glasses  at  1000°C. 
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one  standard  deviation.  Using  the  average  exponential,  an  average  value 
of  the  pre-exponential  for  each  glass  was  calculated.  Utilizing  these 
resultst  an  empirical  equation  for  Cqys|Kt  as  a  function  of  weight  percent 
dissolved  substrate  (P)  and  temperature  was  derived. 

C0r.-|KT  *  f - - -  ]  exp  (-49626/T)  (13) 

0  sl  1  [  1  +  .172  P  +  .00686  Pz  J 

The  Cqysi Ky  values  obtained  from  the  slopes  of  the  experimental  lines  (Figs. 

3-5)  and  from  Eq.  13  are  compared  in  Fig.  6. 

Knowing  CQ  from  the  RUO2  solubility  measurements  [3]  and  taking  ysj 
2 

to  be  20  yJ/cm  ,  Ky  was  calculated  as  a  function  of  composition  and  tempera- 

2 

ture.  The  y$j  value  of  20  pJ/cm  was  suggested  [7]  to  represent  the  solid- 
liquid  interfacial  energy  of  the  Mg0-Ca0-Si02  system.  The  calculated  values 
of  KT,  listed  in  Table  I,  appear  to  go  through  a  maximum  near  4  w/o  substrate 
doping  for  any  given  temperature.  This  anomalous  behavior  is  probably  due 
to  the  assumption  that  the  interfacial  energy  is  independent  of  glass  compo¬ 
sition.  Ru02  solubility  results  [3]  showed  that  the  activity  of  Ru02  in 
glass  increases  rapidly  with  glass  composition  up  to  2  w/o  substrate  doping; 
the  interfacial  energy  may  exhibit  a  composition  dependence  similar  to  that 
of  the  activity  of  Ru02-  Although  the  transfer  coefficient  could  not  be 
determined  precisely  from  the  ripening  results,  the  product  Cqysi Ky  is  suf¬ 
ficient  to  predict  Ru02  sintering  based  on  a  liquid-phase  sintering  model 
presented  in  the  following  section. 
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TABLE  I 


Variation  of  Kj  with  Composition  and  Temperature 


Composition 

Kj  (m/sec)  x  10' 

^ - 

888°C 

950°C 

1000°C 

Standard  Glass 

5.61 

28.3 

191.0 

4  w/o  Sub.  Glass 

6.55 

31.2 

229.0 

10  w/o  Sub.  Glass 

3.45 

16.0 

107.0 

14  w/o  Sub.  Glass 

2.13 

— 

76.0 

4.  Ru02  SINTERING 


Consider  a  pair  of  spherical  particles  in  contact,  completely  immersed 
in  a  liquid.  If  r  is  the  radii  of  the  particles,  the  concentration  of 
species  at  the  solid-liquid  interface  away  from  the  neck  region  (C.j)  rela¬ 
tive  to  that  at  a  flat  interface  (C0)  can  be  written  as: 


In 


Ci  2yslM 
Cq  '  rpRT 


If  C-j  not  much  different  from  CQ,  this  equation  can  be  approximated  as: 

2y 

C1  "  Co  “  rpRT 


'si  M  Co 


04) 


If  x  and  y  are  the  principal  radii  of  curvature  of  the  neck  as  shown  in 
Fig.  7,  the  concentration  of  species  in  the  neck  region  (C2)  relative  to 


Figure  7.  Geometry  of  Neck  Between  Particles  Undergoing 
Liquid-Phase  Sintering. 


A 
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that  at  a  flat  interface  is  given  by 


In 


'2  .  Ysl  * 
'  ’  pRT 


Y$i  M 


for  x  »  y 


Coe- 


rsl 


M  C 


ypRT 


(15) 


In  these  equations,  M  is  the  molecular  weight  of  the  solid  and  p  is  the 
density  of  the  solid.  Combining  Eqs.  14  and  15  gives 


AC  =  C-j 


Ysl  M  Co  T  2 
pRT  L  r 


AC  = 


Ysl  M  Co 
ypRT 


for  r  »  y 


(16) 


If  the  neck  growth  between  the  particles  is  limited  by  the  interface 
reaction,  the  rate  of  mass  transfer,  J(gm/s),  can  be  written  as 


J  =  KtAaC  (17) 

2 

in  which  Kj  (cm/s)  is  the  transfer  coefficient,  A(cm  )  is  the  interfacial 
area,  and  AC(gm/cm  )  is  the  concentration  difference  between  the  free-surface 
and  the  neck  area.  For  transport  limited  by  redeposition  at  the  neck,  the 
interfacial  area.  A,  in  Eq.  17  is  the  neck  area  and  from  Fig.  7, 

A  =  1/2  (2wy )  ( 2irx ) 

=  2n2xy 


Substituting  for  A  and  AC  in  Eq.  17  gives 


By  conservation  of  mass. 


J  =  d  (Vnp)/dt 


where  Vn  is  the  neck  volume.  From  Fig.  7, 


and 


Vn  «  1/2  (2 y) 


y  =  (r  -  a)  =  r(l  -  cos  e)  *  2r  V  “  §7 


2  l 

T  “  2r 


Therefore, 


J  -  d  (^~)  p/d t 


_  2irpx" 
r 


dx/dt 


(19) 


Combining  Eqs.  18  and  19  gives 


„2j..  _  1tKTys1MCo  r  dt 

X  CIX  ~  n 

p% 


which  integrates  to 


#x,3  _  3ltKTYslMCor  2t 

(?)  7m 


(20) 


If  the  neck  growth  is  controlled  by  the  diffusion  of  species  through 
the  liquid,  the  rate  of  mass  transfer  is  given  by 


,  D  A  AC 
J  "  0.78  r 


(21) 


in  which  D(cnr/s)  is  the  diffusion  coefficient  of  the  slowest  moving  species 
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in  the  liquid,  and  0.78  r  is  the  average  diffusion  length  which  is  taken 
to  be  equal  to  1/8  of  the  periphery  of  the  projected  particle.  Substituting 
for  A  and  AC  in  Eq.  21,  and  by  proceeding  as  for  the  interface  controlled 
situation,  it  can  be  shown  that 


The  effects  of  composition  and  temperature  on  the  ripening  of  Ru02  in 
glass  was  given  by  Eq.  13.  By  assuming  that  the  transfer  coefficient  is  the 
same  for  the  ripening  and  the  neck  growth  processes,  the  kinetics  of  neck 
growth  between  Ru02  particles  can  be  computed  by  combining  Eqs.  20  and  13. 
The  resulting  equation  is 


3  Mrt  17.29  x  10  exp  (-49626/T) 
P2RT  (1  +  0.172P  +  0.00686P2) 


The  time  required  for  the  initial  stage  sintering  (x/r  =  0.3)  of  Ru02 

particles  of  different  sizes  in  the  standard  glass  is  given  as  a  function 

0 

of  temperature  in  Table  II.  At  800  C,  for  particles  of  5  nm  (50  A)  radius, 
it  takes  about  165  seconds  to  complete  the  initial  stage  sintering;  whereas 
for  particles  of  1  pm  radius,  it  takes  nearly  2  1/2  months  to  achieve  the 
same  ratio  of  neck  to  particle  radius.  The  effect  of  temperature  on  the 
neck  growth  process,  as  indicated  by  the  high  activation  energy  for  ripening 
is  tremendous.  For  example,  the  time  required  by  5  nm  particles  to  complete 
the  initial  stage  sintering  is  increased  by  a  factor  of  100  if  the  tempera¬ 
ture  is  lowered  from  800  to  700°C. 
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TABLE  II 

Calculated  Initial  Stage  Neck  Growth  Kinetics 
of  RUO2  in  the  Standard  Glass 


Temperature 

Time(s)  for  x/r  =  0.3 

B 

QJQQjjl 

c 

650 

2.60  x  105 

1.04  x  106 

1.04  x  108 

1.04  x  1010 

1.04  x  1012 

700 

1.73  x  104 

6.92  x  104 

-- 

-- 

-- 

750 

1.50  x  103 

6.01  x  103 

— 

— 

— 

800 

1.65  x  102 

6.58  x  102 

-- 

— 

— 

850 

1.65  x  102 

6.58  x  102 

— 

-- 

— 

900 

3.49 

1.39  x  10 

1.39  x  103 

1.39  x  105 
- 1 

1.39  x  107 
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The  effect  of  glass  composition  on  the  neck  growth  kinetics  is  indicated 

by  Eq.  23.  Additions  of  4,  10,  and  14  w/o  substrate  to  the  standard  glass 

increase  the  time  required  for  the  completion  of  initial  stage  sintering  by 

factors  of  approximately  2,  3.5,  and  5,  respectively.  Equation  23  also 

indicates  that  for  a  given  glass  composition,  temperature,  and  time,  the 

neck  size  would  be  proportional  to  the  1/3  power  of  the  particle  size.  The 

neck  size  calculated  as  a  function  of  particle  size  and  sintering 

temperature  for  a  10  minute  sintering  in  the  standard  glass  is  given  in 

Table  III.  A  comparison  of  Tables  II  and  III  indicates  that  although 

larger  particles  take  longer  times  (proportional  to  the  square  of  their 

size)  to  complete  the  initial  stage  sintering  (x/r  =  0.3),  they  can  achieve, 

1  /3 

for  a  given  time,  a  larger  neck  radius  (proportional  to  r  '  )  than  the 
smaller  particles. 
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TABLE  III 


Dependence  of  Neck  Size  on  Particle  Size  and  Temperature 


Temperature 

c 

Neck  Size 

(x)  in  nm  for  t  =  10  min. 

— ■ 

IJQQSES 

n m 

mQQjjj 

650 

0.199 

0.250 

0.539 

1.161 

2.502 

700 

0.490 

0.617 

1.330 

2.865 

6.174 

750 

1.106 

1.394 

3.002 

6.468 

13.94 

800 

>1.5* 

2.914 

6.277 

13.52 

29.14 

850 

>1.5* 

>3.0* 

12.28 

26.46 

57.02 

900 

>1.5* 

>3.0* 

22.68 

48.86 

105.2 

*  Neck  growth  calculations  using  Eq.  23  are  valid  for  x/r  <0.3  only. 
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5.  FUTURE  WORK 

All  experimental  work  on  this  contract  has  been  completed.  The  results 
from  studies  of  glass  viscosity,  surface  tension  and  density,  substrate  dis¬ 
solution  rate,  and  Ru02  solubility  and  ripening  will  be  correlated  utilizing 
the  previously  developed  models  for  microstructure  development,  and  the 
influence  of  glass  composition  established.  The  effects  of  substrate  dis¬ 
solution  on  charge  transport  processes  in  nonsintered  contacts  will  be 
modeled,  and  the  dependence  of  both  the  glass  properties  and  the  electrical 
properties  of  the  nonsintered  contacts  on  glass  composition  will  be  incor¬ 
porated  into  a  revised  charge  transport  model  for  thick  film  resistors. 
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7.  STATEMENT  OF  ESTIMATED  COSTS 


Contract  No.  N0019-79-C-0240 
February  1,  1979  -  January  31,  1980 


$  40,000.00 
32,079.00 

$  7,921.00 

Planned  Expenditures  (Approximate) 

November  $3,420 

December  $2,275 

January  $2,225 


Beginning  Fund  Balance 
Funds  Expended  Through  10/31/79 

Funds  Remaining 


